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We test for the existence of a spin-glass phase transition, the de Almeida-Thouless line, in an externally- 
applied (random) magnetic field by performing Monte Carlo simulations on a power-law diluted one- 
dimensional Ising spin glass for very large system sizes. We find that an Almeida-Thouless line only occurs 
in the mean field regime, which corresponds, for a short-range spin glass, to dimension d larger than 6. 
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Perhaps the most surprising prediction of the mean-field 
theory of spin glasses is that an Ising spin glass has a line 
of transitions in an external magnetic field, the de Almeida- 
Thouless (AT) 111] line. This instability line separates a 
high-temperature high-field paramagnetic phase where relax- 
ation times — ^possibly very large — stay finite, from a low- 
temperature low-field phase where the energy landscape has 
valleys separated by truly infinite barriers in the thermody- 
namic limit. The AT line, an ergodic to non-ergodic transition 
with no change in symmetry, is perhaps the most striking pre- 
diction of the mean-field theory of spin glasses. Whether or 
not it occurs in realistic systems is a major unsolved problem. 

The existence or otherwise absence of an AT line in real 
(short-range) spin glasses is also a key feature distinguishing 
the two most popular scenarios for the nature of the spin- 
glass state below the (zero-field) transition temperature: the 
replica-symmetry breaking (RSB) picture of Parisi 101, and 
the "droplet picture" of Fisher and Huse flQl- The RSB pic- 
ture assumes that the behavior of real spin glasses is very sim- 
ilar to that of the mean-field solution yj] of the Sherrington- 
Kirkpatrick infinite-range model. Since the mean-field model 
has a stable spin-glass state in a field and thus has an AT line, 
it is proposed that this also occurs for any short-range system 
with a finite temperature transition in zero field. By contrast, 
the droplet picture makes certain assumptions about the na- 
ture of the low-energy, large-scale excitations (droplets) from 
which one finds no AT line in any dimension. 

Experimentally, it has been harder to determine if an AT 
line occurs than to show that there is a transition in zero field. 
For the latter case the divergence of the nonlinear susceptibil- 
ity provides a clear signature of the transition. Unfortunately, 
the nonlinear susceptibility does not diverge in a field, i.e., 
along the AT line. However, as noted by two of us iH there is 
a closely-related static quantity which diverges on the AT line 
and which can be measured in simulations, albeit not in ex- 
periments. A finite-size scaling analysis of the two-point cor- 
relation length indicated the absence of an AT line for three- 
dimensional (3D) Ising spin glasses IHI^I- Subsequently, the 
same idea was applied to a one-dimensional (ID) model in 
Ref. 101 (referred to from now on as KY), in which every spin 



interacts with every other spin in the system with a strength 
which falls off with a power of the distance. By varying the 
power, one can simulate the whole range of possible behav- 
iors ||4i|7l|8|], from infinite-range, through mean field, to non- 
mean field and finally to the absence of a finite-temperature 
transition. This is analogous to changing the space dimension 
d of short-range finite-dimensional models. KY found that 
an AT line does occur for parameter values corresponding to 
the mean-field case (for short-range systems that would be for 
d > 6), but not in the non-mean-field case (d < 6). The pos- 
sibility of a critical dimension above which the AT line occurs 
had been considered before, see for example the discussion in 
Ref. H. 

Model and Obserx'ables. — The model studied by KY is 
fully connected so the CPU time for one Monte Carlo sweep 
(MCS) grows as 0{L^), where L is the number of spins. This 
is inefficient for large L. Recently, this difficulty was removed 
in an elegant way in Ref. US by diluting the interactions and 
fixing the connectivity z. We thus study: 
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where Si = ±1 are Ising spins evenly distributed on a ring 
of length L in order to ensure periodic boundary conditions. 
The sum is over all spins on the chain and the couplings Jij 
are normally distributed with zero mean and standard devia- 
tion unity (independent of distance). The dilution matrix 
takes values 1 or 0, and a nonzero entry appears with probabil- 
ity pij, where p^j ^ r^^" with Ty = (-L/tt) sin(7r|i - i\lL) 
representing the geometric distance between the spins. The 
power (7 is a key parameter of the model. To avoid the prob- 
ability of placing a bond being larger than 1, a short-distance 
cutoff is applied and thus we take 



= 1 - exp(-A/r: 
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The constant A is determined numerically by fixing the aver- 
age coordination number z. Note that this model has the same 
long range interactions on average, [J£]av ^ ^l^lj ^ in KY, 
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but has only Lz/2 bonds rather than L{L — l)/2. Hence the 
hnear scaUng of the CPU time for one MCS. 

As in the fully-connected case by varying a one can 
tune the model in Eq. ([TJ from the infinite-range to the short- 
range universality class. For < cr < 1/2 the model is in the 
infinite-range universality class in the sense that the parameter 
A vanishes for iV — > oo, and for tj = it corresponds to the 



scaling. 



Viana-Bray model ll ill . For 1/2 < a < 2/3 the model de- 
scribes a mean-field long-range spin glass, corresponding — in 
the analogy with short-range systems — to a short-range model 
in dimension above the upper critical dimension d > du = Q 
1 1211 . For 2/3 < (T < 1 the model has non-mean-field critical 
behavior with a finite transition temperature Tc- For cr > 1, 
the transition temperature is zero. We are interested in finite- 
range models which have a non-zero Tc, i.e. 1/2 < cr < 1. 

A rough correspondence between a value of cr in the long- 
range ID model and the value of a space dimension in a 
short-range model can be obtained from 
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where ri{d) is the critical exponent jy for the short-range 
model, which is zero in the mean-field regime. Equation (|3]l 
has the following required properties (i) d oo corresponds 
to cr 1/2, (ii) the upper critical dimension du = 6 cor- 
responds to cru = 2/3, and (iii) the lower critical dimension, 
which is where di—2 + ri{di) — 0, corresponds to cr; = 1. For 
example, in 3D, rj ~ 0.384(9) ifisll and thus the corresponding 
exponent is a ~ 0.90. 

In this study we set the average coordination number to 
Zav = 6 and use site-dependent random fields hi chosen from 
a Gaussian distribution with zero mean [/i^Jav = and stan- 
dard deviation [hf]^'^ = H^. The latter has the advantage 
that we can perform a detailed test for equilibration of the 
data when using Gaussian-distributed interactions SQ (see 
below). 

To determine the existence of an AT line, we compute the 
two-point finite-size correlation length Hill El. For this 
we start by determining the wave-vector-dependent spin-glass 
susceptibility given by 
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where (• • • )t denotes a thermal average and [• • • ]av an aver- 
age over the disorder. To avoid bias, each thermal average 
is obtained from a separate copy of the spins, so we simu- 
late four copies at each temperature. The correlation length is 
given by 1711 
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where = 27r/L is the smallest non-zero wave-vector com- 
patible with the boundary conditions. According to finite-size 



^L/L^X[L'/''iT-T,)] , 
a/i"/'^ A'[Li/3(T-r,)] 



cr > 2/3 

1/2 < cr < 2/3, (6) 



with 1/ = l/(2cr — 1) in the mean-field regime ISfl. Hence, if 
there is a transition at T = Tc, data for /L (^l /L"^^ in the 
mean field region) for different system sizes L should cross at 
Tc. 

We also present data for xsG = Xsg(O), which has the 
finite-size scaling form 



XSG T^/^C[Li/3(r - Tc)] , 



(T > 2/3 
1/2 <cr< 2/3. (7) 



Hence curves of xsg/T^ (Xsg/T^/"^ in the mean-field 
regime) should also intersect. This is particularly useful for 
long-range models since rj is given by the naive expression 
2 — ?7 = 2cr — 1 exactly. 

As discussed in KY, for the simulations to be in equilibrium 
with Gaussian fields and bonds, the following equality must 
hold; 
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where q = ^J(S'i)^]av is the spin overlap, qi = 
N^^ X^i j ^ijiSiSj)"^ is the link overlap of a given sample, 
and Nb is the number of nonzero bonds of the sample. To 
speed up equilibration we use theparallel tempering (ex- 
change) Monte Carlo method ll7l llSll . Simulations are per- 
formed at zero field, as well as at = 0.1, a value consid- 
erably smaller than Tc(i?R = 0) for the values of cr studied. 
For details see Table H] 

Results. — We start by showing in Fig. [Ha) data for £,l/L 
against T for cr = 0.75 in zero field, for several system sizes. 
The data intersect cleanly at Tc ~ 1.50 indicating a transition 
at that point, see Eq. (|6]l. The inset shows xsg/T^^'' using 
the exact value rj = 1.5. 

In contrast to Fig.llla), which shows the expected zero-field 
transition for cr = 0.75, Fig.lHb) shows no intersections in a 
small field = 0.1 [approximately 0.067 of the zero-field 
Tc shown in Fig. [Ha)]. Thus there is no AT line for cr = 0.75, 
except possibly for even smaller values of the field. Note that 
cr — 0.75 is in the non-mean-field regime (2/3 < cr < 1). 
Whereas the data for cr — 0.75 for small sizes merge, and it 
is only for the larger sizes that the data do not even meet, for 
cr = 0.85 — deeper in the non-mean-field regime — even the 
data for small sizes do not meet at any temperature down to 
T = 0.30, see Fig. [nc). 

For comparison we also show data in the mean-field regime 
where an AT line is expected to occur For cr = 0.60 and 
= 0.1 there is a clear intersection, see Fig. [Hd). The 
temperature of the intersections is slightly different in the two 
cases, about 1.60 for ^l/T^/'^ and about 1.75 for xsg/L^^'^, 
suggesting finite-size effects, possibly due to long ne gativ e 
tails in the spin overlap distribution; see Fig.|2]and Ref. 11911 . 
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FIG. 1: (Color online) Panel (a) Finite-size correlation length divided by L as a function of T for different sizes for J/r = and o — 0.75 
(non-mean-field region). The inset shows xsg/L^~'' using the exact value ri — 3 — 2a = 1.5. In both cases the data cross indicating a 
phase transition at zero field. Panel (b): Same as (a) but for Hn — 0.1. The absence of an intersection down to low T shows that there is no 
transition in a field [the shaded area corresponds to Tc{Hn = 0)]. The inset shows data for a bimodal (± J) distribution of bonds, as used in 
Ref. flgjl . for sizes L = 256 to 1024 on a linear topology. While Ref. flgfl find a finite-temperature transition (shaded area in the inset) we see 
no sign of it. The absence of a transition is even more clear in panel (c) where we show data as in (b) but for a — 0.85, i.e., deeper into the 
non-mean-field regime. In panel (d) we show data for the correlation length divided by L"/"^ (= L^^^) as a function of T for different sizes 
for Hy{, = 0.1 and ct = 0.60 (in the mean-field region). The inset shows xsg/L^^'^- The intersections show that there is a transition in a field, 
i.e., an AT line for this value of a. 



We note that very recent work by Leuzzi et al. fl^l comes 
to a different conclusion. Using Eq. ([T]) with bimodally- 
distributed disorder they find a transition in a field in the non- 
mean-field regime, in particular for a = 0.75 and H-r, = 0.1, 
where we do not find a transition, see Fig. [Ub). We have 
no explanation for this discrepancy. We have done several 
checks, including developing two versions of the code inde- 
pendently and verifying that they give the same results. Fur- 
thermore, we have simulated the model with the same bimodal 



disorder and geomeh"y as used in Ref. lfl9ll . as well as the same 
field and a values, finding no signature of a transition [see the 
inset toFig.Hb)]. 

Summary. — Our conclusion, based on numerical results, 
is that there is an "upper critical dimension" close to 6 for the 
AT line. This agrees with KY but disagrees with Ref. lfl9ll . 
This conclusion is distinct from RSB theory which predicts 
an AT line in any space dimension with a zero-field transition, 
and the droplet picture flSl^ according to which there is no 
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no. 2: (Color online) Distribution of the spin overlap q for a — 
0.75, T = 0.71 and Hn = 0.1. Even for the largest L studied there 
is a tail which extends into the negative-g region. 

TABLE I: Parameters of the simulations for different field strengths 
Hr and exponents a. N^a. is the number of samples, A^sw is the 
total number of Monte Carlo sweeps, Tmin is the lowest temperature 
simulated, and Nt is the number of temperatures used in the parallel 
tempering method for each system size L. The last column shows 
the parameter A [Eq. l|2j] fixing Zav = 6 neighbors. 



a 




L 






T ■ 

^ mm 


Nt 


A 


0.60 


0.10 


128 


8000 


8192 


0.480 


46 


0.99458 


0.60 


0.10 


256 


8000 


32768 


0.480 


46 


0.90363 


0.60 


0.10 


512 


5000 


131072 


0.480 


46 


0.83827 


0.60 


0.10 


1024 


5000 


524288 


0.480 


46 


0.78926 


0.60 


0.10 


2048 


4500 


65536 


1.393 


26 


0.75140 


0.75 


0.00 


128 


5000 


32768 


0.300 


50 


1.71141 


0.75 


0.00 


256 


5000 


32768 


0.300 


50 


1.64289 


0.75 


0.00 


512 


5000 


524288 


0.300 


50 


1.59859 


0.75 


0.00 


1024 


2900 


2097152 


0.300 


50 


1.56903 


0.75 


0.00 


2048 


1000 


2097152 


0.480 


46 


1.54892 


0.75 


0.00 


4096 


1000 


65536 


1.192 


31 


1.53506 


0.75 


0.00 


8192 


500 


131072 


1.192 


31 


1.52544 


0.75 


0.10 


128 


5000 


32768 


0.480 


46 


1.71141 


0.75 


0.10 


256 


5000 


131072 


0.480 


46 


1.64289 


0.75 


0.10 


512 


5000 


262144 


0.480 


46 


1.59859 


0.75 


0.10 


1024 


5000 


524288 


0.480 


46 


1.56903 


0.75 


0.10 


2048 


2800 


524288 


0.710 


39 


1.54892 


0.85 


0.10 


128 


6000 


16384 


0.300 


50 


2.39485 


0.85 


0.10 


256 


6000 


65536 


0.300 


50 


2.34867 


0.85 


0.10 


512 


6800 


524288 


0.300 


50 


2.32189 


0.85 


0.10 


1024 


2500 


2097152 


0.300 


50 


2.30592 



AT line in any finite dimension. Of course the numerical data 
cannot rule out a transition at extremely small fields. 

Note added in proof: We have recently heard (G. Parisi, 
private communication) that there is an error in the analysis of 
Ref. d, and that their results for a = 0.75 are now much 
more similar to ours. 
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